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Abstract 

Background: Extensive and intensive studies on tine unfolding of proteins require appropriate theoretical model and 
parameter to clearly illustrate the feature and characteristic of the unfolding system. Over the past several decades, four 
approaches have been proposed to describe the interaction between proteins and denaturants, but some ambiguity and 
deviations usually occur in the explanation of the experimental data. 

Methodology/Principal Findings: In this work, a theoretical model was presented to show the dependency of the residual 
activity ratio of the proteins on the molar denaturant concentration. Through the characteristic unfolding parameters /c, and 
Am, in this model, the distribution, transition and thermodynamic stability of protein conformations during the unfolding 
process can be quantitatively described. This model was tested with the two-state unfolding of bovine heart cytochrome c 
and the three-state unfolding of hen egg white lysozyme induced by both guanidine hydrochloride and urea, the four-state 
unfolding of bovine carbonic anhydrase b induced by guanidine hydrochloride and the unfolding of some other proteins 
induced by denaturants. The results illustrated that this model could be used accurately to reveal the distribution and 
transition of protein conformations in the presence of different concentrations of denaturants and to evaluate the unfolding 
tendency and thermodynamic stability of different conformations. In most denaturant-induced unfolding of proteins, the 
unfolding became increasingly hard in next transition step and the proteins became more unstable as they attained next 
successive stable conformation. 

Conclusions/Significance: This work presents a useful method for people to study the unfolding of proteins and may be 
used to describe the unfolding and refolding of other biopolymers induced by denaturants, inducers, efc. 
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Introduction 

Studies on the unfolding and refolding of proteins have been 
attracting the attention of researchers in the fields of biochemistry, 
biophysics and structural biology for the past decade, and they 
play an important role in gaining an understanding of the 
relationship between the structure of proteins and their functions 
[1-5]. To date, the conformational transition of proteins mainly 
originates from studies on the unfolding and refolding of proteins 
induced by denaturants, pH, heat, etc [6,7]. Of these factors, 
denaturants may be the most important. Through the conforma- 
tional transition of the proteins induced by denaturants, valuable 
information has been obtained on the self-organization, function, 
misfolding and aggregation of proteins [8-10]. 

In fact, the denaturant-induced conformational transition of 
proteins can be considered the result of the interaction between 
the protein and denaturant molecules. Through studies on 
conformational transitions, it is now believed that globular 
proteins are usually unfolded through one or several denatured 



but compact equilibrium intermediate(s) [11,12], and they may 
exist in four different types of conformation states: native states (N), 
molten globule states (MG), pre-molten globule states (PMG), and 
completely unfolded states (U) [13,14]. One of the most 
characteristic features of the partially folded conformations (MG 
and PMG) is their combination of properties that are typical of the 
N and U states. Proteins in MG states have a globular structure 
that is typical of native globular proteins, and they may have a 
native-like secondary structure and a native-like folding pattern; 
while proteins in PMG states have no rigid tertiary structure and 
are characterized by a considerable secondary structure, although 
it is much less pronounced than that of N or MG state proteins. 
Proteins in this state are considerably less compact than those in 
the MG or N states, but they are still more compact than those in 
the random coil [15-19]. 

Over the past several decades, four approaches have been 
proposed to describe the interaction between proteins and 
denaturants. All of these approaches use the following general 
relationship to describe the unfolding free energy (AGli) of the 
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protein: 

AG = AG°^-AGex (a) 

where AGu" indicates the unfolding free energy of the protein in 
the absence of denaturants and AGex denotes the difference in the 
excess free energ\' of the interaction between the denaturant and 
the native or denatured protein in its conformational transition. 
The first approach is the purely phenomenological linear 
extrapolation method shown in Eqn. (b): 

AGu = AGO -m-[D] (b) 

where m indicates the slope of the dependence of AGU on the 
molar denaturant concentration ([D]) in denaturation solution. 
Therefore, the unfolding free energy (AGy") of the protein in the 
absence of denaturants can be obtained by the linear extrapolation 
of the AGu values observed in the denaturation transition zone 
when the denaturant concentration (|D]) is zero [20]. The second 
approach relies on the free energy of transferring amino acids and 
peptides from water into the aqueous denaturant solution [21-23]. 
These data can be used in conjunction with the amino acid 
composition of the protein to derive a measure of the change in the 
solvent-accessible surface area with differing denaturant concen- 
trations [23]. The third approach views the denaturation of the 
protein as the result of the binding of the denaturant by the native 
and denatured protein [23]: 

AGu = AC^-AnWTlnil+kd) (c) 

where k indicates the binding constant of the protein to the 
denaturant, An denotes the difference in the number of binding 
sites on the surface of the native and denatured protein, and R and 
T represent the gas constant and the absolute temperature, 
respectively. The fourth approach is similar to the third approach. 
In this method, a binding formalism is used in its mathematical 
development, but it treats the act of denaturation more like 
solvation than binding [24—27]. A key distinction between the 
solvent exchange approach and the denaturant binding approach 
is the formal inclusion of the solvent in the interaction between the 
protein and denaturant in the former. Of these four approaches, 
the most often used is the linear extrapolation approach. 

However, it should be realized that although these approaches 
have been used to describe the conformational transition of 
proteins induced by denaturants, some ambiguity and deviations 
usually occur in the explanation of the expc-rimcntal data. For 
example, for the linear extrapolation method, although this 
method has been widely used to derive the unfolding free energy 
(AGu") of the protein in the absence of denaturants [28-32], the 
following problems usually exist, (i) The method is semi-empirical 
and lacks an authentic theoretical basis. The linear relationship 
between the denaturant concentration ([D]) and the unfolding free 
energy (AGu) of the protein is only an experimental phenomenon; 
it is not based on any authentic theory and is not certain to exist 
over the whole denaturant concentration range. Therefore, 
whether the linear extrapolation method is reasonable to describe 
the interactions between proteins and denaturants is a controver- 
sial issue, (ii) Theoretically, the conformational transition of 
protein should be described by the thermodynamic equilibrium 
constant of its conformational transitions rather than by the 
concentration ratio {C\/ Cg) of its two adjacent conformation states 
(A and B). For a given unfolding system, when temperature and 
denaturant concentration are known, the unfolding free energy 



(AGu) of the proteins only depends upon the stability of its two 
adjacent conformation states and has no relation to the denaturant 
concentration. However, the concentration ratio {C\/C^) in the 
presence of different denaturant concentrations is used in this 
method and it is involved not only in the temperature and the 
characteristics of the denaturant but also in the denaturant 
concentration, (iii) Some inconsistent and even paradoxical 
conclusions frequendy occur using this semi-empirical method. 
For example, in most cases, the unfolding free energy (AGl") of the 
protein in the absence of denaturant, which is derived from one 
kind of denaturant, is very different from another type of 
denaturant. Similarly, some ambiguity and deviations also exist 
in the other three methods. For example, the values of parameter 
m derived from the denaturant binding method are always 
approximately 10-100% less than those derived from calorimetric 
studies of urea and guanidine hydrochloride solutions [33]. 
Therefore, it is necessary and urgent to develop appropriate 
theoretical model and characteristic parameter to accurately 
describe the interaction between the protein and denaturant and 
to clearly show the features of the unfolding system. 

In this work, a theoretical model is presented to illustrate the 
distribution, transition and thermodynamic stability of protein 
conformations during the denaturant-induced unfolding. This 
model shows the relationship between the residual activity ratio of 
the proteins and the denaturant concentration. Using this model, 
two characteristic unfolding parameters, and Am;, which 
separately indicate the thermodynamic equilibrium constant for 
the unfolding of the proteins and the change in the number of the 
denaturants associated with each protein molecule between stable 
conformation states, can be simultaneously derived. Furthermore, 
using these two characteristic unfolding parameters, the distribu- 
tion and transition of each stable conformation of the protein can 
be well described and the unfolding tendency and thermodynamic 
stability of each conformation can be clearly predicted over the 
unfolding processes. 

Theory 

Interactions between proteins and denaturants. For the 

sake of convenience, we assumed that (i) as shown in Figure 1, 
denaturant-induced protein unfolding usually is a progressive 
procedure. The protein is generally unfolded from its native state 
(Njj) to its completely unfolded state (Ub) through one or more 
intermediate states (I_d,) in a denaturation solution. Although the 
unfolding of the protein as a whole may occur through non- 
equilibrium thermodynamics over a wide denaturant concentra- 
tion range, the protein may be in local thermodynamic 
equilibrium between adjacent conformational transitions under a 
local denaturant concentration area. Therefore, the association- 
dissociation interaction between the protein and denaturant can 
still be designated as thermodynamic equilibrium, (ii) An 
interaction always exists between the protein and denaturants 
during unfolding. The given protein may be in the state, lo 
state or U/; state in the denaturation solution, and when the 
protein exists in different stable conformation states it can associate 
with different numbers of denaturant molecules (Figure 2). 
Therefore, when they are unfolded from one stable conformation 
state to another adjacent state, a change must also occur in the 
number of the denaturant molecules associated with the protein. 

Therefore, the progressive unfolding of the protein from Nj) 
state to Ui3 state through n Id, state can be expressed as: 
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Figure 1. Denaturant-induced unfolding of a protein from its native state to completely unfolded state through n intermediate 
states. Nd, Id, and Uq denote the native state, intermediate state and completely unfolded state of a protein molecule, respectively, k, indicates the 
thermodynamic equilibrium constant for the unfolding of the protein from one stable conformation state to the next. The symbol " < ' " implies 
that the unfolding of the protein by the denaturant is only at a local thermodynamic equilibrium under a local denaturant concentration range. 
doi:1 0.1 371/journal.pone.0091 1 29.g001 
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where D denotes the denaturant; ^ (i= 1, 2, n+1) represents 
the thermodynamic equilibrium constants for the unfolding of the 
protein between stable conformation states; and Am; (i= 1, 2, 
n+1) denotes the change in the number of denaturant molecules 
that are associated with each protein molecule between stable 
conformation states. 
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In Eqn. (1), A:; (i= 1, 2, n+1) can be expressed as: 



Ami ' fx 1 vj-w/^ni') 



[Iz)„_,l-[D1 



[Uz)] 



(2) 



+ 1 



where [N^j], [Id-] (i= 1, 2, n) and [Uj)] indicate the 
equilibrium concentrations (mol/L) of the native state, interme- 
diate state and completely unfolded state of the protein in the 
denaturation solution, respectively, and |D] denotes the denatur- 
ant concentration (mol/L). 
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0 0 0 0 0 0 
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(aJ and (aq) : a protein molecule associated without or with a number of 
denaturant molecules under a stable conformation state A; 

expand CS>: a protein molecule associated without or with a number of 
denaturant molecules under a stable conformation state B; 
Q : a denaturant molecule in denaturation solution. 

Figure 2. Interactions between the protein and denaturant molecules in the denaturant-induced unfolding of proteins. 

doi:1 0.1 371 /journal.pone.0091 1 29.g002 
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From Eqn. (2), we can obtain in sequence: 



(3) 



[lD„]=k„ilo„_,]m^"=kyk2 fc„.[D]^l+^"2 + - + A-«^ 

[U„]=^„+r[Iz,„l-[Dl^'"«+i 

= kvk2 A:„+r[D]^i+^"'2+ 



If we define tlie residual activity ratio (r) of tlie protein at a given 
denaturant concentration as tlie ratio of [N^i] to their initial total 
concentration (Cj) in the denaturation solution, then we have: 



[Nz,] 



+ ] + [!/)„] + [Uz,] 



(4) 



In Eqn. (4), C, is the sum of [N/,], [I/jJ and [U/J in the 
denaturation solution. Therefore, substituting [I^),] (i= 1, 2, n) 
and [UiJ from Eqn. (3) in Eqn. (4) and then re-arranging it, we 
have: 



l+A:r[D]*»'i +k,-k2iD\^'' + ■ ■ • +ki-k2- ■ ■ ■ •/c„+i-[D]*"i+^2+ 



(5) 



Eqn. (5) shows the dependency of the residual activity ratio of 
the protein on the denaturant concentration. 

Derivation of characteristic unfolding parameters ki and 
Ami, For the sake of mathematical processing convenience, we 
wiU use the three-state protein unfolding mechanism as an 
example to show the derivation of the characteristic unfolding 
parameters k; and Am-,. In this case, Eqn. (5) can be simplified to: 



1 



(6) 



where ki and k2 indicate the thermodynamic equilibrium constants 
for the unfolding of the protein from N^, to I^j and from Ij^ioVj), 
respectively. Ami ^md Amg denote the changes in the number of 
the denaturant molecules associated with a protein molecule 
during the unfolding from N^j to Ir, and from Ij^ to Vr,, 
respectively. 

In the denaturant-induced three-state unfolding of proteins, it is 
not possible that each of stable conformation state is evenly 
distributed at different denaturant concentrations. However, it is 
possible that the protein mainly exists in one or two stable 
conformation states in a certain denaturant concentration. 
Therefore, Eqn. (6) can be flirther simplified. 

(i) When the denaturant concentration is relatively low, the 
protein mainly exists in N^, state and Ijj state. Thus, IJ £, state 
can be omitted. In this case, Eqn. (6) can be simplified to: 



1 



(7) 



when the residual activity ratio is not equal to zero or 1 and 
the denaturant concentration is not equal to zero. Re- 
arranging Eqn. (7) and taking the natural logarithm on both 
sides of the equation, we have:. 



ln( - - 1) = Inki + Ami -Inp] 



(8) 



As shown in Eqn. (8), in this situation, a linear relationship exists 
between ln(l/r— 1) and ln[D]. Therefore, by determining the 
residual activity ratio over a range of relatively low denaturant 
concentrations, we can derive the two characteristic unfolding 
parameters ki and Ami from the intercept and slope of the linear 
regression plots of ln(l/?— 1) versus ln[D]. 

(ii) When the denaturant concentration is relatively high, the 
protein mainly exists in In state and \Jj) state. Thus, N/) state 
can be ignored. In this situation, Eqn. (6) can be simplified to: 



1 



A:r[D]^i(l-hA:2-[D]'^2) 



(9) 



when both the residual activity ratio and the denaturant 

concentration are not equal to zero. Re-arranging Eqn. (9) 
and taking the natural logarithm on both sides of the 
equation, we have: 



ln( 



1 

r-kv[U]' 



-1)= lnik2 + Aw2Tn[D] 



(10) 



As shown in Eqn. (10), a linear relationship also exists between 



ln[l/(r,i;i-[D]'^ 



1] and ln|T)] in this case. Therefore, by 



determining the residual activity ratio over a range of relatively 
high denaturant concentrations, we can further derive the two 
characteristic unfolding parameters k2 and Am^ from the intercept 



and slope of the linear regression plots of ln[l/(r^i-[D] 



■1] 



\Trsus In[D]. 

In conclusion, by determining the residual activity ratio of the 
protein over a range of denaturant concentrations, we can 
separately derive the characteristic unfolding parameters ki, k^, 
Anil and Am2 in the three-state unfolding of proteins induced by 
denaturants. 

Distribution and transition of protein conformations. In 

this section, we wiU still use the three-state unfolding model as an 
example. In the denaturant-induced three-state unfolding of 
proteins, at a given denaturant concentration, the molar 

fractions/Nfl , /iD^nd /uoOf the native state, intermediate state 
and completely unfolded state can be respectively expressed as: 



[Nz)l 



M = 

Q [Nz,] + [Ifl] + [Uz,]' 

f _ [M _ M 

Ct [Nz,] + [Iz,] + [Uz,]' 



(11) 



ho'- 



C, [Nz,] + [Iz,] + [Uz,] 



[Uz,] 



PLOS ONE I www.plosone.org 



4 



March 2014 | Volume 9 | Issue 3 | e91129 



Transition and Stability of Protein Conformations 



Obviously, the following relationship exists among the molar 
fractions/Nj, , /ij, and/uj, : 



(12) 



By substituting [I/)] and \Ud\ from Eqn. (3) in Eqn. (1 1) and 
then re-arranging Eqn. (11), we have: 



ln(^ -l) = lnA:i+AOTi ln[D], 

( , r^A^n. - 1 ) = In ^2 + Am2 • In [D] , 



(18) 



ln( 



1)= lnA;3-l-Am3-ln[D] 



fin = 



1 



ytr[D]^i 



(13) 



l+fcr[D]'^i+A:i-A:2-[D]'^i+'^"'2 



Eqn. (13) shows the effect of the denaturant concentration on 
the molar fractions. When the denaturant concentration ap- 
proaches zero,_/Nj,is close to l,_/ij,and/uj,are close to zero: 



lim /Nn = l> lim /in = 1™ /un=0 (14) 

[01-0 -° [D]^0 ^ [D]^0 ^ ^ 



When the denaturant concentration approaches infinity, fijc is 
close to l,fif„&ndfij, are close to zero: 



According to Eqn. (13), as long as the characteristic unfolding 
parameters ki, ^2, Ami 'ind Am2 'ire given, we can derive the molar 
fractions /no, /icand/u^ at different denaturant concentrations, 
and furthermore, we can describe the distribution and transition of 
these diflFerent conformation states. 

Similarly, for the denaturant-induced two-state unfolding of 
proteins, the following hnear relationship should exist between 
ln(l/r-l) and ln[D]: 



ln(^-l)=lnA;-hAwln[D] 



(16) 



Through the linear regression plot of ln(l/?— 1) versus ln|D], 
the two characteristic unfolding parameters k and Am can be 
derived from the intercept and slope, respectively. In addition, at 
different denaturant concentrations, the molar fractions fti^zxid 
/uj, can be derived from the following equation: 



fvo = 



[Nz)] 
Ct 



1 



[Nz,H-[Uz,l 



1+^[D]"' 

,1 Am 



(17) 



Ct [Ni,] + [Uz,] l+;t-[D]' 



For the denaturant-induced four-state unfolding of proteins, the 
characteristic unfolding parameters ki, k2, kj, Ami, Am^ and Ams 
can also be derived using the following hnear regression equations: 



At different denaturant concentrations, the molar fractions 
_/nj,i/iipj i/iipjand /u^ can be derived from the following equations: 



M = 

C, [Nz,]-F[Ifli]-l-[Ic2]-F[Uc] 
1 

l-HA:r[D]^'i +/crfc-[D]^"'i+^'"2+fc,.A;2-/t3-[D]'^i+'^2+'^3 ' 



no 



1 C, [Ni,]-h[Ifli]-l-[Ifl2] + [Uz)] 



/li). 



(19) 



'2 Ct [Nfl]-h[Ii>j]-K[Ii,2] + [Ui,] 

;ti-/t2-[D]^i+^"2 



l-t-ztrfDl^i -FA:r/t2-[Dl^i+^2+A:,.A;2-A:3-[Dl^i+'^2+Am3 ' 



[Uz.] 



M = 

C, [Nc]-h[Iz),] + [lD2] + [Ui,] 

' \\kv\Of"^ +A:r/t2-[Dl^l+'^2+;tr/k2-A:3-[Dl'^i+'^2+Am3 



Results 

Two-State Unfolding of Bovine Heart Cytochrome c 

As shown in our previous work, the unfolding of the bovine 
heart cytochrome c induced by both guanidine hydrochloride and 
urea is a typical two-state process [34]. Besides, in our previous 
work, the structure changes of bovine heart cytochrome c during 
the two unfolding processes were also studied. From the results, it 
can be found that in the absence of denaturants, about 27% of the 
Trp residues in bovine heart cytochrome c are exposed to the 
quencher molecules; when the urea concentrations are about 1.25 
and 3.0 mol/L, about 54% and 97% of the Trp residues are 
separately exposed to the quencher; while when the guanidine 
hydrochloride concentrations are about 3.0 and 6.0 mol/L, about 
51% and 99% of the Trp residues are separately exposed. All the 
evidence indicated that the tertiary structure of bovine heart 
cytochrome c is destroyed with the increase of denaturants 
concentration, the molecule structure become loose and all the 
tryptophan residues originally embedded in the interior of their 
molecules are exposed to the surface of their molecules. 
Meanwhile, the size-exclusion chromatographers and native 
polyacrylamide gel electrophoreses results showed that during 
the denaturants-induced unfolding, bovine heart cytochrome c 
molecules exist only in a unimolecular form and their poly- 
molecular aggregates are not formed all along. 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-1 0123456789 

Denaturant concentration/mol L"^ 

Figure 3. Residual activity ratios [i) of bovine heart cytochrome c exposed to different concentrations of guanidine hydrochloride 
or urea. □: guanidine hydrochloride; A: urea. The concentration of bovine heart cytochrome c was 0.50 mg/mL, and the experimental temperature 

was 25"C. 

doi:1 0.1 371/journal.pone.0091 1 29.g003 
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Figure 4. Plots of ln(1/r— 1) vs. ln[D] for the unfolding of bovine heart cytochrome c induced by guanidine hydrochloride and urea. 

O: guanidine hydrochloride; □: urea. 
doi:1 0.1 371/journal.pone.0091 1 29.g004 
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Table 1. Regression correlation coefficients (R^) and characteristic unfolding parameters k and Am for the unfolding of bovine 
heart cytochrome c Induced by guanldine hydrochloride and urea. 



denaturant concentration 

Denaturants conformational transitions Ayl mol^ Am range 

Guanidine hydrochloride ., Jj^,, 0.998 2.41 x10"' 4.85 0-6.0 L mol"' 

Urea ., , , 0.998 3.79x10"^ 5.53 0-8.0 L mol"' 



doi:l 0.1 371 /journal.pone.0091 1 29.t001 

In the solution containing different concentrations of denatur- 
ants, tlie residual activity ratio of bovine heart cytochrome c was 
determined (Figure 3). For the two unfolding processes of bovine 
heart cytochrome c, the residual activity ratio continuously 
decrease from 100% to nearly 0% as the urea concentration 
gradually increase from 0.0 to approximately 6.0 mol/L or the 
guanidine hydrochloride concentration gradually increase from 
0.0 to approximately 3.0 mol/L. Except when the urea concen- 
tration is equal to either 0.0 or 6.0 mol/L, the residual activity 
ratio of protein induced by urea is higher than that induced by 
guanidine hydrochloride. Therefore, it can be inferred that bovine 
heart cytochrome c molecules are more easily unfolded with 
guanidine hydrochloride than with urea. 

For the two unfolding processes of bovine heart cytochrome c, 
linear regression plots of ln(l / r— 1) vs. ln|D] were created (Figure 4) 
according to Eqn. (16) and the linear regression correlation 
coefficients {R ) and the characteristic unfolding parameters k and 
Am were derived (Table 1). From these calculations, we 
determined that the good linear relationship exists between ln(l/ 
r—l) and ln|D] for both the two unfolding processes and that the 
regression correlation coefficients {R^) are not lower than 0.998. 

As shown in Figures 5A and 5B, based on the above 
characteristic unfolding parameters k and Am and by using Eqn. 
(17), the denaturant-induced distribution and transition of native 
and completely unfolded bovine heart cytochrome c can be 
visualized by plotting their molar fractions /iM^and fij^as the 
function of the guanidine hydrochloride or urea concentration 
[D] . These figures clearly show the denaturant-induced unfolding 



landscape of bovine heart cytochrome c, the percentages of N/) 
state and Vo states in the presence of different concentrations of 
denaturants, and how state gradually transforms into U^i state 
during the unfolding processes. For example, from Figures 5A and 
5B, we know that (i) both guanidine hydrochloride- and urea- 
induced unfolding of bovine heart cytochrome c are two-state 
processes, (ii) In the urea-induced unfolding of bovine heart 
cytochrome c, when no urea molecules exist in the solution, all of 
the protein exists in its native state. When the urea concentration is 
1.0 mol/L, approximately 99% of the protein stiU exists in N^j 
state, and only approximately 1% exists in \Jj) state. When the 
urea concentration reaches approximately 5.7 mol/L, only 1% of 
the protein remains in state, and 99% exists in U/) state. When 
the urea concentration reaches approximately 6.0 mol/L, all of 
the bovine heart cytochrome c exists in state, (iii) In the 
guanidine hydrochloride-induced unfolding of bovine heart 
cytochrome c, when the guanidine hydrochloride concentration 
is approximately 0.20 mol/L, only 1% of the protein exists in "Nq 
state, and 99% exists in U/j state. When the guanidine 
hydrochloride concentration is approximately 1.5 mol/L, approx- 
imately 50% of the protein is in N/) state, and 50% is in Uj) state. 
When the guanidine hydrochloride concentration reaches approx- 
imately 3.0 mol/L, all of the protein exists in state. 

Three-State Unfolding of Hen Egg White Lysozyme 

The unfolding of hen egg white lysozyme induced by both 
guanidine hydrochloride and urea has been shown to be a three- 
state process [36]. In previous work, the results of intrinsic 




Guanidine hydrochloride concentration/mol L" 



Urea concentration/mol L" 



Figure 5. Molar fractions /nz, and /y^ of native and completely unfolded bovine heart cytochrome c exposed to different 
concentrations of guanidine hydrochloride (A) or urea (B). O: native bovine heart cytochrome c; □: completely unfolded bovine heart 
cytochrome c. 

doi:1 0.1 371/journal.pone.0091 1 29.g005 
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Figure 6. Residual activity ratios [i) of hen egg white lysozyme exposed to different concentrations of guanldlne hydrochloride or 
urea. A: guanidine hydrochloride; C: urea. The concentration of hen egg white lysozyme was 0.50 mg/mL, and the experimental temperature was 
25°C. 

doi:1 0.1 371/journal.pone.0091 1 29.g006 



fluorescence anisotropy of lysozyme induced by denaturants 
indicated that, by raising the denaturants concentration, the 
flexibility and freedom degree of the six Trp residues in lysozyme 
molecules increase, which meant that the micro-environment 
where Trp residues locate in transform to a state with more loose, 
more flexible and smaller space occupy structure. This work also 
demonstrated that the deactivation precede the structure change 
of lysozyme molecule. 

In the solution containing different concentrations of denatur- 
ants, the residual activity ratio of hen egg white lysozyme was 
determined (Figure 6). For the two unfolding processes of hen egg 
white lysozyme, the residual activity ratio continuously decrease 
from 100% to almost 0% as the urea concentration gradually 



increase from 0.0 to approximately 8.0 mol/L or the guanidine 
hydrochloride concentration gradually increase from 0.0 to 
approximately 6.0 mol/L. Therefore, hen egg \vliit(; lysozyme is 
also more easily unfolded with guanidine hydrochloride than with 
urea. 

For the two unfolding processes of the hen egg white lysozyme, 
the linear regression plots of ln(l/>— 1) vs. ln[D] and ln[l/ 
(rAi-[D]'^i)-l] vs. ln[D] were created (Figures 7A and 7B) 
according to Eqns. (8) and (10) and their linear regression 
correlation coefficients (R^) and the characteristic unfolding 
parameters ki, k^. Ami and Am2 were derived (Table 2). From 
these calculations, we determined that a good linear relation- 
ship exists between ln(l/r— 1) and ln[D] and between 




Figure 7. Plots of ln(1//--1) vs. ln[D] (A) and ln[1/(A''Ari '[D]'^'">)-1] vs. ln[D] (B) for the unfolding of hen egg white lysozyme Induced 
by guanidine hydrochloride and urea. A: guanidine hydrochloride; O: urea. 
doi:l 0.1 371/journal.pone.0091 1 29.g007 
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Table 2. Regression correlation coefficients (R^) and characteristic unfolding parameters /c,, /c2, Am, and for the unfolding of 
hen egg white lysozymes induced by guanidine hydrochloride and urea. 



denaturant concentration 

Denaturants conformational transitions Ai/L mol^ Am, range 

Guanidine hydrochloride ., _hMn^ 0.986 4.16x10"^ 3.41 0-2.5 L mol"' 
' Np < > Id 

, 1 1 0.989 3.64x10"' 3.86 3.5-6.0 L mol"' 

Urea ., JnAm 0.985 1.78x10"' 5.14 0-2.8 L mor' 

No < > Id 

, 1 1 0.990 2.95x10"^ 2.58 4.5-8.0 L mol"' 
Id < ' Uq 



doi:1 0.1 371 /journal.pone.0091 1 29.t002 

ln[l/(rA:i-[D]^"")-l] and lii[D] and that all of their regression 
correlation coefficients are not lower than 0.985 for this process. 

The molar fractions , /ij, and /uj, of hen egg white lysozyme 
at different concentrations of denaturants were derived (Figures 8 A 
and 8B) according to Eqn. (13). These results indicate that (i) both 
the unfolding of guanidine hydrochloride- and urea-induced hen 
egg white lysozyme are three-state processes, and the protein 
molecule can only be transformed from its state to Vo state 
through an intermediate state, (ii) In the guanidine hydrochloride 
solution, the hen egg white lysozyme molecule is first unfolded 
from its N^, state to In state as the guanidine hydrochloride 
concentration increases to approximately 3.0 mol/L and then is 
further unfolded from I^, state to Vo state as the guanidine 
hydrochloride concentration further increases to approximately 
6.0 mol/L. The maximum molar fraction of the protein in If, state 
is approximately 0.42. (iii) In the urea solution, the egg white 
lysozyme molecule is first unfolded from its state to In state 
when the urea concentration increases to approximately 4.0 mol/ 
L and then is further unfolded from In state to U^j state when the 
urea concentration further increases to approximately 8.0 mol/L. 
The maximum molar fraction of hen egg white lysozyme in In 
state is approximately 0.40. 



Four-State Unfolding of Bovine Carbonic Anhydrase b 

The guanidine hydrochloride-induced unfolding of bovine 
carbonic anhydrase b has been shown to be a four-state process 
[38]. To the structure change of carbonic anhydrase molecule 
during the guanidine hydrochloride-induced unfolding, it was 
found that the first intermediate has no rigid tertiary structure but 
is almost as compact as the native protein, and it meets all the 
usual requirements of the molten globule state. The second 
intermediate, which is less compact than the molten globule state 
but much more compact than the unfolded state, represents a 
novel, "pre-molten globule" state. The CD results indicated that 
the secondary structure of carbonic anhydrase is destroyed in two 
phases: first, upon transition between molten globule and pre- 
molten globule states and then upon a complete unfolding of 
protein molecules, besides, both the two intermediates have 
substantial secondary structure. The ANS binding results suggest- 
ed that the pre-molten globule state binds ANS about five times 
more weakly than the molten globule state. And the size-exclusion 
chromatography results illustrated that the pre-molten globule 
state is reasonably compact. 

The residual activity ratio of bovine carbonic anhydrase b 
exposed to different concentrations of denaturants was determined 
now (Figure 9). For the two unfolding processes of bovine carbonic 
anhydrase b, the residual activity ratio of protein continuously 
decrease from 100% to nearly 0% as the urea concentration 




1 2 3 4 6 

Guanidine hydrochloride concentration/mol L"' 



Urea concentration/mol L" 



Figure 8. Molar fractions f^^, fi^ and f\j„ of native, intermediate and completely unfolded hen egg white lysozyme exposed to 
different concentrations of guanidine hydrochloride (A) or urea (B). □: native egg white lysozyme; O: intermediate egg white lysozyme; A: 
completely unfolded egg white lysozyme. 
doi:1 0.1 371/journal.pone.0091 1 29.g008 
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■\ I 1 I I I I I 1 I 1 I I 
0 1 2 3 4 5 

Denaturant concentration/mol L"^ 

Figure 9. Residual activity ratios (i) of bovine carbonic anhydrase b exposed to different concentrations of guanidine 
hydrochloride or urea. O: guanidine hydrochloride; □: urea. The concentration of bovine carbonic anhydrase b was 0.50 mg/mL, and the 
experimental temperature was 25'C. 
doi:10.1371/journal.pone.0091129.g009 



gradually increase to approximately 4.4 moI/L or the guanidine 
hydrochloride concentration gradually increase to approximately 
2.2 mol/L. Therefore, bovine carbonic anhydrase b is also more 
easily unfolded with guanidine hydrochloride than with urea. 

Similarly, for the guanidine hydrochloride-induced four-state 
unfolding of bovine carbonic anhydrase h, plots of the residual 
activity ratio of protein vs. the denaturant concentration were also 
created according to Eqn. (18), and their linear regression 
correlation coefficients [R ) and the characteristic unfolding 
parameters ki, k2, k^, Ami, ^'"2 "ind Am^ were also derived 
(Table 3). As shown in Table 3, all of the regression correlation 
coefficients are not lower than 0.985. 

Based on Eqn. (19), we can also derive the molar fractions /n^ , 
/mGc/pmGz) and /u^ of bovine carbonic anhydrase b molecules at 
different concentrations of guanidine hydrochloride (Figure 10). 
Through these calculations, we found that (i) the unfolding of 
guanidine hydrochloride-induced bovine carbonic anhydrase /> is a 
typical four-state process and that bovine carbonic anhydrase b 



can only be transformed from its state to Vj) state through two 
sequential intermediate states, the Img state and the I/a/f; state, (ii) 
In the guanidine hydrochloride solution, bovine carbonic anhy- 
drase b is first unfolded from its N^j state to MGf, state as the 
guanidine hydrochloride concentration increases to approximately 
1.4 mol/L. The protein is then further unfolded from MG/j state 
to its PMGx) state as the guanidine hydrochloride concentration 
further increases to approximately 1.7 mol/L. Finally, bovine 
carbonic anhydrase b transforms from its FMGd state to Uij state 
as the guanidine hydrochloride concentration further increases to 
approximately 2.2 mol/L. The concentrations at which the 
number of molecules in the molten MG^, state and the PMG^) 
state reaches its maximum are approximately 1.40 and 1.7 mol/L 
guanidine hydrochloride, respectively. These results are approx- 
imately consistent with those from a previously published study 
[38]. 



Table 3. Regression correlation coefficients (R^) and characteristic unfolding parameters k^, l<2, k^, Ami, and Am^ for the 
unfolding of bovine carbonic anhydrase b induced by guanidine hydrochloride. 



guanidine hydrochloride 

Conformational transitions AiA mor^ \m, concentration range 

ti.Aim 0.988 1.35x10"^ 10.02 0-1.3 L-mor' 

fa An; 0.997 6.91 x10"" 14.70 1.4-1.8 Lmor' 

MGd ( > PMGd 

pi^Q fa''^") y 0.992 7.71 x10"^ 17.30 1.9-2.4 L mor' 



doi:1 0.1 371 /journal.pone.0091 1 29.t003 
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Figure 10. Molar fractions /ns'/mgz,i/pmgz, and /u„ of native, molten globule, pre-molten globule and completely unfolded bovine 
carbonic anhydrase b exposed to different concentrations of guanidine hydrochloride. O: native bovine carbonic anhydrase b; □: 
molten globule bovine carbonic anhydrase b; A: pre-molten globule bovine carbonic anhydrase b; O: completely unfolded bovine carbonic 
anhydrase fa. 

dol:10.1371/journal.pone.0091129.g010 



Unfolding of other Proteins 

In all of the unfolding processes described above, we were able 
to derive the characteristic unfolding parameters ki and Ami by 
exploiting the relationship between the residual activity ratio (r) 
and the denaturant concentration ([D]). Furthermore, through 
these characteristic unfolding parameters, we were able to describe 
the distribution and transition profiles of each stable conforma- 
tional state at different denaturant concentrations. However, we 
can also derive the characteristic unfolding parameters and Am; 
through the distribution and transition data of the denaturant- 
induced unfolding of the proteins. 

For example, in the denaturant-induced three-state unfolding of 
hen egg white lysozyme, when the denaturant concentration is 
relatively low, the protein mainly exists in its N/) state and !/> state, 
and \Jb state can be omitted. In this case, using Eqn. (13), we 
obtain: 



ln(- l)=lnA:i-l-Awiln[D] 



(21) 



When the denaturant concentration is relatively high, hen egg 
white lysozyme mainly exists in its In state and U/) state, and its N^j 
state can be ignored. In this case, through Eqn. (13), we obtain: 



fx 



1 



(22) 



Similarly, when the molar fraction /u^ of egg white lysozyme in 
state is not equal to zero or 1 and by re-arranging Eqn. (22) 
and taking the natural logarithm, we also obtain: 



1 



(20) 



ln( - 1) = \nk2 + Aw2 ■ln[D] 
Jvd 



(23) 



When the molar fraction /jsf^ of egg white lysozyme in Njj state 
is not equal to zero or 1 and by re-arranging Eqn. (20) and taking 
the natural logarithm, we have: 



Eqns. (21) and (23) show that in the denaturant-induced three- 
state unfolding of hen egg white lysozyme, a linear relationship 
exists between ln(7^ 1) and ln[D] and between ln(7^ 1) and 

ln[D]. From their intercepts and slopes, the characteristic 
unfolding parameters ki. Ami and k^, Am2 can be simultaneously 
derived. 
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Table 4. Regression correlation coefficients (R^) and characteristic unfolding parameters /c, and Am, for the unfolding of proteins 
induced with different denaturants. 



Conformational transitions 



Ai/L mol ' 



A/77i 



denaturant concentration range 



The three-state unfolding of /Mactoglobulin induced by methanol 

0.990 



.An 17 

MGn / > Un 



0.988 



8.61 xlO 
1.51 xlO" 



9.17 
5.39 



The four-state unfolding of the apo form of myristoylated NCS-1 induced by guanidine hydrochloride 



At. Awn 



0.997 
0.995 



3.50x10 
1.78x10" 



9.39 
9.65 



The four-state unfolding of the apo form of non-myristoylated NCS-1 induced by guanidine hydrochloride 



r>l2, 



0.998 



The four-state unfolding of adenylated tslig induced by guanidine hydrochloride 

0.996 4.92x10"^ 



kj .Ann 



PMGn-Z=>Un 



0.991 
0.987 



4.29x10 
9.92x10" 



The four-state unfolding of de-adenylated tslig induced by guanidine hydrochloride 

0.988 1.13x10"' 



ki.Ami 
Nn ( > MG, 



MGd ( " > PIVIGd 



Pi\/lGo-zz=>UD 



0.991 
0.994 



3.06x10 
6.68x10" 



The four-state unfolding of /f-lactamase induced by guanidine hydrochloride 

0.987 1.99x10' 



k\ .Ami 
Np ( > MG, 



MGn 



I>PIVlGn 



k-f ,A/«i 
PMGo , > Un 



0.987 
0.989 



5.29x10' 
4.25 



4.72 
9.29 
13.62 

3.53 

16.79 

13.74 

2.66 
2.17 
4.70 



The quasi-five-state unfolding of rabbit muscle creatine kinase induced by guanidine hydrochloride 

0.986 2.77x10^ 4.70 



f ''^"'> MGo 

ko.Ami 

MGd ( > PMGd 

kj ,A/«3 

PMGp -p >Hc 



0.991 
0.988 
0.995 



7.95x10^ 

5.30x10"' 

2.53x10"" 



8.99 
9.69 
12.91 



0-25% (V/V) 
50-70% (VA/) 





0.991 


5.55 


2.52 


0-1.5 L-mol 











1.5-4.0 L mol 
4.5-7.0 L mol 



ki.Ami 


0.993 


1.44x10' 


4.50 


0-1.0 L-mol 











1.0-3.5 L mol 



l2. 


0.997 


4.57x10"^ 


10.40 


4.8-7.0 L mol 











0-2.5 L mol 

2.5- 3.4 L mol 

3.6- 4.6 L mol 

0-2.5 L mol 
2.5-4.0 L mol 
5.0-7.0 L mol 

0-0.3 L mol 
0.5-0.8 L mol 
1.0-2.0 L mol 

0-0.2 L-mol"' 
0.3-0.6 L-mol"' 
0.6-2.2 L-mol"' 
2.2-5.8 L-mol"' 



doi:1 0.1 371/journal.pone.0091 1 29.t004 



Using the distribution and transition data for the unfolding of 
other protein molecules induced by denaturants in the literature, 
as shown in Table 4, we derived their regression correlation 
coefficients {R^) and the characteristic unfolding parameters and 
Am;. These data included the three-state unfolding of [i- 
lactoglobulin induced by 0-80% (v/v) methanol [41], the four- 
state unfolding of the apo form of myristoylated NCS-1 induced by 
0-6.0 mol/L guanidine hydrochloride and the apo form of non- 
myristoylated NCS-1 induced by 0-5.0 mol/L guanidine hydro- 
chloride [42] , the four-state unfolding of adenylated tslig induced 
by 0-5.0 mol/L guanidine hydrochloride and de-adenylated tshg 
induced by 0-7.0 mol/L guanidine hydrochloride [43], the four- 
state unfolding of /^-lactamase induced by 0-2.0 mol/L guanidine 
hydrochloride [44] and the quasi-five-state unfolding of rabbit 
muscle creatine kinase induced by 0-6.0 mol/L guanidine 



hydrochloride [45] . As shown in Table 4, all of their regression 
correlation coefficients are greater than 0.985. 

Discussion 

Fitting of the Theoretical Model 

In all of the above fitting analyses, including that of Eqn. ( 1 6) to 
the two-state unfolding data of bovine heart cytochrome c, Eqns. 
(8) and (10) to the three-state unfolding data of hen egg white 
lysozyme, Eqn. (18) to the four-state unfolding data of bovine 
carbonic anhydrase b and Eqns. (21) and (23) to the three-state 
unfolding data of hen egg white lysozyme, at least 10 sets of 
experimental data were used to fit these equations and aU of the 
regression correlation coefficients were greater than 0.985. 
Therefore, it can be concluded that this theoretical model can 
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Figure 11. Plots of InA; vs. Ami for the unfolding of some proteins induced by denaturants. O: guanidine hydrochloride; □: urea; A: 
methanol. 

doi:10.1371/journal.pone.0091 129.g01 1 



be used to fit the residual activity ratio and the distribution and 
transition data of proteins exposed to difierent denaturant 
concentrations with a 99% degree of confidence. 

Characteristic Unfolding Parameters k| and Ami 

As shown in the mathematical expression of the theoretical 
model, Eqn. (5) includes two characteristic unfolding parameters 
and Am;. Both and Am; only depend on the inherent 
characteristics of the protein and denaturant molecules but not 
relate to the denaturant concentration, k, is the thermodynamic 
equilibrium constant for the denaturant-induced unfolding of 
proteins, and it shows the unfolding tendency of protein from one 
stable conformation state to the next. Am; is the difference in the 
number of denaturant molecules associated with a protein 
molecule between stable conformation states, and it shows the 
ability of a denaturant to unfold a protein in a given unfolding 
system. Therefore, it can be inferred that a relatively larger k^or a 
relatively smaller Am; means that it is easier for the protein 
transition from one stable conformation state to the next; whereas 
a relatively smaller ^ or a relatively larger Am; means that it is 
harder for this transition to occur. 

For example, as shown in Table 1 for the two-state unfolding of 
bovine heart cytochrome c, the characteristic unfolding parameter 
k is 3.79x10^'' and 2.41 xlO~' L-moF' for urea and guanidine 
hydrochloride, respectively. The value for urea is approximately 
64 times that for guanidine hydrochloride. Meanwhile,the Am is 
5.53 and 4.85 for urea and guanidine hydrochloride respectively. 
Both k and Am show that the guanidine hydrochloride-induced 
conformational transition of bovine heart cytochrome c from its 



native state to its completely unfolded state is easier than the urea- 
induced transition. 

In another example, as shown in Table 2 for the three-state 
unfolding of hen egg white lysozyme induced by guanidine 
hydrochloride, the characteristic unfolding parameters ki and ^2 
are 4. 16x10 and 3.64x10 L-mol , respectively. The value 
of ki is approximately 1 1 times that of ^2- Meanwhile, Am; and 
Am2 are 3.41 and 3.86 respectively. These results show that the 
conformational transition of hen egg white lysozyme induced by 
guanidine hydrochloride from its native conformation state (N/j) to 
its intermediate conformation state (I^)) is easier than from its 
intermediate conformation state (In) to its completely unfolded 
conformation state (Ub). For the three-state unfolding of hen egg 
white lysozyme induced by urea, the characteristic unfolding 
parameters ki and k2 are 1.78x10 ' and 2.9,5x10 L-mol 
while Ami and Am2 are 5.14 and 2.78, respectively. These results 
show that the urea-induced conformational transition of hen egg 
white lysozyme from N/j to In is harder than the transition from 
to Vo. However, for the three-state unfolding of the hen egg white 
lysozymes, the total characteristic unfolding parameter k^, which is 
equal to the product of the characteristic unfolding parameters ki 
and ^2) is 1.51 xlO~* and 5.25x10"^ L-mol~' for guanidine 
hydrochloride and urea, respectively. While the total characteristic 
unfolding parameter Amt, which is equal to the sum of the 
characteristic unfolding parameters Am; and Am2, is 7.27 and 7.72 
for the two unfolding processes, respectively. Both of these 
parameters show that the conformational transition of hen egg 
white lysozyme from its native stable conformation state to its 
completely unfolded state induced by guanidine hydrochloride is 
easier than that induced by urea. 
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And as shown in Tables 3 and 4, similar situations also exist for 
the unfolding of other proteins induced by denaturants. 

In addition, through carefully analyzing and comparing ^ and 
Am; in Tables 1 to 4, we found that a reverse change exists 
between the characteristic unfolding parameters in aU of these 
unfolding systems, in other words, in each step of the- conformation 
state transition, a relatively large is always followed by a 
relatively small Am;, and vice versa. From the above results, we can 
further predict whether a linear regression relationship definitely 
exists between the characteristic unfolding parameters and Am;. 
However, as shown in Figure 1 1 and Eqn. (24), when linear 
regression plots of In^ vs. Am; were created, no statistically linear 
relationship existed between the characteristic unfolding param- 
eters: 

\nki = - 1 .49Ami + 1M,R^ = 0.528 (24) 

This result means that k-^ and Am; are two independent variables. 
They separately show the characteristic features of the denaturant- 
induced unfolding of proteins from two different aspects. 

From both the fitting analyses of the mathematical expressions 
to the experimental data and the physical meanings of the 
characteristic unfolding parameters, it may be concluded that this 
theoretical model can be used to effectively describe the unfolding 
of proteins induced by denaturants. 

Unfolding Tendency and Thermodynamic Stability 

The characteristic unfolding parameters ^ and Am; in this 
model are the thermodynamic equilibrium constant for the 
denaturant-induced unfolding of proteins and the change in the 
number of the denaturant molecules that are associated with a 
protein molecule between two successive stable conformation 
states, respectively. Therefore, for a given unfolding system, we 
can estimate the unfolding tendency and thermodynamic stability 
of each of its conformation states. 

In Tables 2, 3 and 4, through carefully comparing the 
characteristic unfolding parameters, it can be found that, in most 
of the examined unfolding processes, the ki in the subsequent 
transition step is smaller than that in the previous step, and the An;; 
is larger in the uc'xt transition stej) than in the pre\'ious transition 
step. However, these trends are not the case for the urea-induced 
three-state unfolding of hen egg white lysozyme; the methanol- 
induced three-state unfolding of ;6-LG; the guanidine hydrochlo- 
ride-induced four-state unfolding of de-adenylated Tslig; and the 
guanidine hydrochloride-induced four-state unfolding of /?-lacta- 
mase. These data indicate that in most denaturant-induced protein 
unfolding processc-s, th(- unfolding of the protein becomes 
increasingly harder with each transition step. 

At a given temperature (T), the free energy change (AGu) in the 
denaturant-induced unfolding of a protein from its stable 
conformation state (A) to its adjacent stable conformation state 
(B) can be expressed as: 

AGu=-RT'lnA:, (25) 

where R indicates the gas constant, and ki denotes the 
thermodynamic equilibrium constant for the unfolding from its 
state A to its adjacent state B. Using Eqn. (25), we can estimate the 
stability of the protein in each stable conformation states during 
the unfolding process. If k^ is more than 1, then the free energy 
change AGu is negative, indicating that the unfolding of the 
protein from state A to its adjacent state B is spontaneous and that 



the protein is more stable in state B than in state A. If A; is less than 
1, then the free energy change AGu is positive, indicating that the 
unfolding of the protein from state A to state B is not spontaneous 
and that the protein is more stable in state A than in state B. If k, is 
close to 1, then the free energy change AGu is close to zero, 
indicating that the unfolding of the protein from state A to state B 
is approximately equal to its refolding from state B to state A and 
that the state A is approximately as stable as state B. 

As shown in Tables 1 to 4, in most denaturant-induced 
unfolding processes, k; is less than 1, and consequentiy, AGu is 
positive. This is the case ("xcept for the guanidine hydrochloride- 
induced four-state transitions of the apo-NCS-1 and four-state 
transitions of the apo-non-myristoylated NCS-l(N/)->Iij3); the 
guanidine hydrochloride-induced four-state transition of )S- 
lactamase(N/)->MG/)->PMGi3-»Uz)); and the guanidine hydro- 
chloride-induced quasi-five-state transition of rabbit muscle 
creatine kinase(Njj->N^->MG/), N^:dimeric state with a rigid 
tertiary structure). Therefore, in most denaturant-induced protein 
unfolding processes, the protein is more unstable in its next 
adjacent stable conformation state than in its previous stable 
conformation state. 

In summary, a theoretical model was presented to show the 
dependence of the residual activity ratio of proteins on the molar 
denaturant concentration in the denaturant-induced unfolding, 
the distribution, transition, unfolding tendency and thermody- 
namic stability of protein conformations can be quantitatively 
described by the two characteristic unfolding parameters ki and 
Am; in this model. In most denaturant-induced unfolding processes 
of proteins, the unfolding becomes increasingly harder in next 
transition step and the protein becomes more unstable as it attains 
each successive stable conformation. This work presents a useful 
method for people to study the unfolding of proteins and may be 
used to describe the unfolding and refolding of other biopolymers 
induced by denaturants, inducers, etc. 

Materials and Methods 

Reagents and Chemicals 

Bovine heart cytochrome c, bovine carbonic anhydrase b and 
dried Micrococcus ly.mdeikticus cells were purchased from Sigma. Hen 
egg white lysozyme (20000 units/mg), yeast cytochrome c, reduced 
and oxidized glutathione (highly pure), urea (highly pure) and 
guanidine hydrochloride (highly pure) were purchased from 
Shanghai Sangon Biological and Technological Service Co., Ltd. 
AU other reagents were of analytical grade without further 
purification, and all solutions were prepared in bi-distiUed water 
and passed through a 0.22 Jim filter. 

Unfolding of Bovine Heart Cytochrome c 

In previous work, we determined the intrinsic fluorescence 
emission spectrum, fluorescence phase diagram, and fluorescence 
quenching and deactivation profile for the unfolding of bovine 
heart cytochrome c induced by urea and guanidine hydrochloride 
and found that the unfolding of bovine heart cytochrome c 
induced by both guanidine hydrochloride and urea was a typical 
two-state process [34]. 

According to Darley-Usmar et al. [35], the activity of bovine 
heart cytochrome c at different concentrations of guanidine 
hydrochloride or urea was determined using the Clark oxygen- 
electrode method, in which the oxidase activity was measured by 
determining the oxygen concentration in solution. Yeast cyto- 
chrome c was used as reference to measure the difference in 
voltage (AmV). The characteristic unfolding parameters k and Am 
for bovine heart cytochrome c under different concentrations of 
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guanidine hydrochloride and urea were derived using Eqn. (16), 
and the molar fractions /n^ and_/u^ of the native state (N/,) and 
completely unfolded state (Vj)j were obtained with Eqn. (17). 

Unfolding of Hen Egg White Lysozyme 

Using the "phase diagram" method of fluorescence, Yang et al. 
found that the unfolding of hen egg white lysozyme induced by 
both guanidine hydrochloride and urea was a typical three-state 

process [36]. 

In the study t)f Rozema and GeUman [37], the activity of hen 
egg white lysozyme exposed to different concentrations of 
guanidine hydrochloride and urea was assayed at 25°C by 
calculating the decrease in the absorbance of the reaction solution 
at 450 nm. The reaction was carried out by adding 0. 1 mL of the 
enzyme solution to 1.0 mL of a 0.25 mg/mL Micrococcus 
fysodeikticus suspension in 0.06 mol/L potassium phosphate 
(pH 6.2). One activity unit corresponded to an absorbance 
decrease of 0.0026 per minute. Using different concentrations of 
guanidine hydrochloride and urea, the characteristic unfolding 
parameters ki and Ami, and k2 and Am2 for hen egg white 
lysozyme were derived using Eqns. (8) and (10). In addition, the 
molar fractions /nb, flo and /u^, of the native state (N/)), 
intermediate state (Ijyj and completely unfolded state (U^j) were 
obtained using Eqn. (13). 

Unfolding of Bovine Carbonic Anhydrase b 

Using far UV circular dichroism spectroscopy, fluorescence 
spectroscopy and size-exclusion chromatography, Uversky and 
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